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Abstract 



We reformulate the Gross-Pitaevskii equation with an external parabolic potential as a discrete 
dynamical system, by using the basis of Hermite functions. We consider small amplitude stationary 
solutions with a single node, called dark solitons, and examine their existence and linear stability. 
Furthermore, we prove the persistence of a periodic motion in a neighborhood of such solutions. 
Our results are corroborated by numerical computations elucidating the existence, linear stability 
and dynamics of the relevant solutions. 

1 Introduction 

We address the Gross-Pitaevskii (GP) equation with an external parabolic potential 



where U{X,T) : M x i-^ C is decaying to zero as |x| — > cxd, e G M is the strength of the external 
potential and o" = 1 (a = — 1) is normalized for the defocusing (focusing) cubic nonlinearity. This 
equation is of particular interest in the context of Bose-Einstein condensates, i.e., dilute alkali vapors 
at near-zero temperatures, where dynamics of localized dips in the ground state trapped by the mag- 
netically induced, confining potential V{X) = e^X'^ is studied in many recent papers, see review in 
A question of particular interest concerns whether the localized density dips oscillate periodically 
near the center point X = of the potential V{X). If the motion of a localized dip is truly periodic, 
the frequency of periodic oscillations is to be found [5j, while if the periodic oscillations are destroyed 
due to emission of radiation, the gradual change in the amplitude of oscillations is to be followed for 
sufficiently small e [16j. Numerical simulations show radiation and amplitude changes if the confining 
parabolic potential is perturbed by a periodic potential while no radiation and time-periodic oscillations 
in the case of purely parabolic confinement |18j . 

If cr = 1, a localized dip on the ground state of the GP equation (jl.ip in the formal limit e ^ 
represents the so-called dark soliton of the defocusing nonlinear Schrodinger (NLS) equation, which is 
the reason why we use the term "dark soliton" for a localized solution. Persistence and stability of 
a dark soliton of the defocusing NLS equation in the presence of an exponentially decaying potential 
V{X) was studied in our previous paper [T7], where methods of Lyapunov-Schmidt reductions, Evans 
functions and the stability theory in Pontryagin space were employed. These methods can not be 
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applied to the potential V{X) = e^X"^ since the potential deforms drastically the spectrum of the 
linearized problem: the continuous spectral band at e = becomes an infinite sequence of isolated 
eigenvalues for e 7^ 0. Therefore, we do not use here the limit e — > 0. Moreover, we transform the GP 
equation (jl.ip to the e-independent form 

iut = -]^Uxx + ]jX^u + a\u\^u, (1.2) 

by the scaling transformation x = \X, t = A^T, and u{x,t) = X~^U{X,T) with A = 

Substitution u{x, t) = e~2*~'^^^(j){x) reduces equation (II. 2p to the second-order non-autonomous ODE 

- i</)"(x) + ix2</,(x) + a</.3(x) = (/" + ^i^)^ (1-3) 

where : R M. A strong solution of the ODE (11.30 is said to be a dark soliton if (f){x) is odd on 
E M, has no zeros on 2; E M4., and decays to zero sufficiently fast as |x| — > 00. A classification of 
all localized solutions of the second-order ODE (II. 3p and their construction with a rigorous shooting 
method is suggested in recent work [3]. 

Substitution of m(x, iw{x, t)] reduces equation (jl.2p to the PDE system 

vt = C-W + 2a(j){x)vw + (7{v'^ + w'^)w, (14) 
-wt = C+v + a(f>{x){3v'^ + w"^) + a{v'^ + w^)v, ^ ' ' 

where {v,w) : M x M? and C± are self-adjoint Schrodinger operators in L^(]R) 



= -^dl + ^x'^-^-fi + 3a(l}'^{x), 
jC_ = -^^dl + lx"^ -\- iJL + a(tP{x). 



(1.5) 



Solutions of the PDE ()1.2p are considered in space 

Hi(M) = {n G //^(M) : G ^^(M)} (1.6) 

equipped with the norm 

Hull?,, = / (|n'(x)|2 + (x2 + l)Kx)p)dx. (1.7) 
Similarly, the domain of operators C± in (jl.Sp is defined in space 

W2(K) = {n G //^(M) : x^ueL^iM.)]. (1.8) 

The PDE system (|1.4p is a Hamiltonian system with the standard symplectic structure and the 
Hamiltonian function in the form 

H =]^{v, C+v) + ^ {w, C^w) + (j{(t)v, + w"^) + '^{v'^ + w^, + w"^), (1.9) 

where (•,•) denotes a standard inner product in L^(M). The Hamiltonian function H is bounded if 
v^w ^ 'Hi(M) and is constant in time t. Due to the gauge invariance of the PDE (jl.2p . there exists an 
additional quantity 

Q = 2{(P, v) + {v, v) + {w, w), (1.10) 

which is constant in time t. Global existence of solutions of the initial-value problem associated with 
the PDE (jl.2p in space u £ TCi(M) for all t G M+ has been proved (see Proposition 2.2 in [6j). 
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Considering the linear part of the PDE system (|1.4|) . one can separate the variables in the form 
v{x,t) = f(x)e^*, w{x,t) = w{x)e'^* and obtain the linear problem 

£+f = —Xw, C-w = Xv (1-11) 

for the spectral parameter A G C and the eigenvector {v,w) £ L^(M, C^). Fix /i G M such that a 
stationary solution of the ODE (|1.3p exists and (p € 'Hi(M). Then, the linear problem (jl.lip admits an 
exact solution 

X = ±i : V = (j)'{x), w = ^ix(j){x), (1-12) 

and {v,w) S C^). Additionally, for any /i E M, for which the stationary solution (j){x) is smooth 

with respect to parameter fi, the linear problem (jl.lip admits another exact solution for zero eigenvalue 
A = of geometric multiplicity one and algebraic multiplicity two: 

C^(t>{x) = 0, = -(t>{x). (1.13) 



The main part of this work is devoted to the study of periodic solutions of the PDE system (II. 4|) for 
values of /i near // = 1. This special value corresponds to the second eigenvalue of the linear Schrodinger 
operator 

£ = 48| + i-^-i (1.14) 

with the eigenfunction (/>(x) = exe~^ 1"^ . Here the parameter e is arbitrary in the linear problem and 
it parameterizes the corresponding family of stationary solutions (//, ^{x)) of the nonlinear ODE (jl.3p . 
which bifurcates from the small-amplitude eigenmode [3]. 

A periodic solution of the PDE system (jl.4p bifurcates from the linear eigenmodes v{x) = Se^'^ (j)' {x) 
and w{x) = ^i6e^'^ x(j){x) corresponding to the eigenvalue pair A = iti of the linear problem (jl.lip . Here 
5 and r are two real-valued parameters, which are arbitrary in the linear problem and parameterize the 
corresponding family of periodic solutions {v,w) of the PDE system (jl.4p . An additional parameter a 
comes from the projection of the solution (v, w) to the geometric kernel of the linear problem (jl.lip 
with the eigenmode v{x) = and w{x) = a(p{x). Using this construction, the main result of our paper 
is described by the following theorem. 



Theorem 1 Let e and 6 be sufficiently small and let a, r be arbitrary. There exists a unique family of 
solutions of the ODE ( fj.gj) such that 



exe 



-x2/2| 



Til < Cie 



/327r 



< C2e\ 



(1.15) 



for some e-independent constants Ci,C2 > 0. There exists a family of time-periodic space-localized 
solutions of the PDE system (fi.^j ) such that {v,w) G Hi(M,M^) for all t G M, 



V { x,t -\- 



27r 



v{x, t), w I x,t + 



2ir 



w{x, t), y{x, t) G 



with the bounds 



If (•, t) — 5(j)'{x) cos(ilt -|- r)| 



Hi 



\w{-,t) - 5[x4>{x)sui{Qt + t) + a(j){x)]\\^^ < CieS"^ 



(1.16) 



(1.17) 
(1.18) 



and 1$! — 1| < C^e 5 for some (e, 5) -independent constants 0-^,0^,0^ > 0. 
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The periodic solution of Theorem [T] has four free parameters (e, 6, r, a) which are associated with 
projections to the four eigenmodes (jl.l2p and (jl.lSp of the hnear problem (jl.lip . Parameters r and 
a can be set to zero due to two obvious symmetries of the PDE (\1.2\) : the gauge transformation 
u{x,t) n(x,t)e*", Va G M and the reversibility transformation u{x,t) i— > u{x, —t), Vt G M. 

Although the eigenmodes ()1.12p and (jl.l3p persist for all /i € M, existence of periodic orbits of 
the GP equation (|1.2p is only proved near fi = 1. This is due to the fact that the non-resonance 
conditions n ^ ImAm, Vn, m S N are proved to be satisfied only in this domain, where Am denote other 
isolated eigenvalues of the linear problem (jl.lip with ImA^ > which are different from X = i. Thus, 
resonances do not occur near the value fi = 1. The construction of the periodic orbit is complicated 
due to the existence of translational eigenmodes associated with the double zero eigenvalue A = of 
the linear problem (jl.lip . 

Our main result is in agreement with Theorem 2.1 in [9], where the Newton's law of particle dynamics 
is obtained in a more general context of multi-dimensional confining potentials and general nonlinear 
functions of the GP equation iip = —V'^ip + V{x)ip — f{ip). The Newton's law is derived for parameters 
{a,p) of the solitary wave solution of the unperturbed equation with V{x) = and it takes the form 

d = 2p, p = -W{a). (1.19) 

Adopting our notations for the time variable and the potential function of the GP equation ()1.2p . we 
rewrite the Newton's law (jl.lOp in the explicit form d + a = 0, which recovers the frequency 17 = 1 of 
the periodic solution of Theorem [T] in the linear approximation 5 — > 0. 

There are several differences between results of Theorem 2.1 in ^ and our Theorem [TJ First, the 
Newton's law (|1.19p is valid on finite time intervals and in the limit when the localization length of 
the stationary solution (j){x) is much smaller than the confinement length of the potential V{x). This 
situation corresponds to the original GP equation (jl.ip in the limit e — > 0. Second, the exact periodicity 
is not guaranteed by the periodic solutions of the Newton's law (ll.lOp because of the remainder terms. 
Lastly, the frequency $7 = 1 of the Newton's law is independent of the nonlinear function f{tp) and the 
nonlinear corrections in {a,p). In our case, the result of Theorem [T] is valid for all time intervals, the 
exact periodicity is guaranteed, and the frequency O changes with parameters 5. On the other hand, 
our results are valid in the limit — > 1, which is far from the limit e — > of the GP equation (jl.ip . 

Note that the oscillations of the dark solitons in the GP equation (11. 2p with the frequency = 1 were 
predicted from the Ehrenfest Theorem in much earlier works (see references in [5] and [9j). However, 
it was argued that this frequency is not observed in numerical simulations of the original GP equation 
()l.ip with (7 = 1 for sufficiently small e [U [TBI I18j . It was suggested in these works (see review in 
pTj ) that dark solitons oscillate with a smaller frequency ^ = We will show that both frequencies 
occur in the spectrum of the linear problem (|l.lip in the corresponding limit but the non-resonance 
conditions are not satisfied for either frequency in this limit. 

Our strategy for the proof of Theorem [1] is to use a complete set of Hermite functions and to 
reformulate the evolution problem for the PDE (II. 2p as an infinite-dimensional discrete dynamical 
system for coefficients of the decomposition (Section 2). Existence of stationary solutions 4>{x) of the 
ODE ()1.3p and spectral stability of stationary solutions in the linear problem (jl.lip are studied in the 
framework of the discrete dynamical system (Section 3). The proof of existence of periodic solutions of 
the PDE system (jl.4p relies on construction of periodic orbits in the discrete dynamical system (Section 
4). The analytical results are verified with numerical approximations of solutions of the ODE ()1.3p . 
eigenvalues of the linear problem (jl.lip and solutions of the GP equation (jl.2p (Section 5). Distribution 
of eigenvalues of the linear problem (jl.lip in the limit /x ^ oo for o" = 1 is also analyzed with formal 
asymptotic methods (Appendix A). 
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2 Formalism of the discrete dynamical system 



The set of Hermite functions is defined by the standard expressions ^ : 

Mx) = ^^=^Hn{x)e-^'/\ Vn = 0,1,2,3,..., (2.1) 
y^2"n!0r 

where Hn{x) denote the Hermite polynomials, e.g. Hq = 1, Hi = 2x, H2 = 4rE^ — 2, = 8x"^ — 12x, 
etc. Since the Hermite functions are eigenfunctions of the linear Schrodinger equation 

-lKix) + lx^Mx) = (^n + ^^Mx), Vn = 0,l,2,3,..., (2.2) 

the Sturm-Liouville theory implies that the set of Hermite functions {(/)„(x)}^q forms an orthogonal 
basis in (M) . The normalization coefficients in the expressions (j2.ip ensure that the Hermite functions 
have unit L^-norm, such that 

{(t>n,4>m) = Sn,m, Vn, 771 = 0, 1, 2, 3, .... (2.3) 

We represent a solution u{x,t) of the GP equation ()1.2p by the series of eigenfunctions 

oo 

u{x, t) = e"5* ^ an{t)Mx) (2.4) 

n=0 

where the components (oq, ai, 02, ...) form a vector a on N. When the series representation (j2.4p is 
substituted to the GP equation (jl.2p . the PDE problem is converted to the discrete dynamical system 



icin — na^i (J ^ ^ -f^n, ni, 712, naf^rai 0^712 '^rta ; — 0,1,2,3,..., (2-5) 

(711,712, "a) 

where Kn,ni,n2,nz = {4'n-,(l>ni<t>n2'Pn3)- We shall use a convention to avoid specifying the range of non- 
negative integers (n.i,n2,n3) and n in the summation signs of the dynamical system ()2.5p . Let /^(N) 
be a weighted discrete /^-space equipped with the standard norm 

oo 

||a||| = ^(1 + n)2^|a„p < oo, Vs G M. (2.6) 

n=0 

Since the set {(/)„(x)}^q forms an orthonormal basis in L^(R), we note the isometry ||m||^2 = 11^11^2) so 
that u G L^(R) if and only if a G /^(N). On the other hand, we need an equivalence between the space 
7^1 (M) for the function u{x) and the space Z^(N) for the vector a. In addition, we need to determine 
the domain and range of the vector field of the discrete dynamical system (j2.5p . These results are 
described in Lemmas [1] and [2j 

Lemma 1 Let u{x) = ^ an4>n{x). Then u £ 'Hi(M) if and only i/a G l^^i^)- 

Proof. It follows directly that 

||n|||^^ = / {\u'ix)\^ + {x^ + l)\u{x)\^)dx 



5 



oo oo „ 

^ ^ aniOnj / [(/>^^(x)0'„2(x) + {x^ + l)0„i (x)(/)„2 (x)] dx 



ni=0 n2=0 
oo oo 

= 2 ^ ^ a„ian2(l + n2)((/)ni,(/'n2) 
ni=0 n2=0 
oo 

= 2V(l + n)|a„|2 = 2||a||22 , 
^ — ' '1/2 

n=0 

where the orthogonality relations (|2.3p have been used. □ 

Remark 1 By the same method, one can prove that u G W2(K) if and only if a G (N). 

Lemma 2 T/ie vector field of the dynamical system 112. 5]) maps /^^2(^) ^'^ ^-1/2 (^)' 

Proof. The vector field of the dynamical system (|2.5p is decomposed into the linear f (a) and nonlinear 
o"g(a) parts, where 

fn — IT'O'm 9n — ^ ^ -^n,ni,n2,?i3'^ni^n2^?i3 i — 0,1,2,3,... 

The linear unbounded part satisfies the estimate 

oo 

l|f(a)ll| =E(1 + ^)''"'I««I'^ II^IlL' (2.7) 

n=0 

such that f : /s+i(N) ^ /^(N) for ah s G M. If a G l\i2{^), then s = -\. The nonhnear vector part 
satisfies the estimate 

oo 

n=0 {n\,n2,'nz) (m\,m2,'nnz) 

oo / oo \ 

= j;(l + n)2^|(0„n,|n|2)|'< ^(l + n)2«||tx0„||22 ||^z||i4 

n=0 \ra=0 / 

< (f;(i+n)2^ii0„iiij Nii^, 

\n=0 / 

oo 

where n(x) = ^ an<pn{x) and all (/>n(a;) are real-valued. By the main theorem of [7], there exists a 

n=0 

constant C > such that 



h l|4 ^ ^log(l + n 

vi + 



<C ^ ; S Vn = 0,1,2,... (2.^ 



Therefore, the series Y1'^=q{^ + "')^*ll</'nlli4 converges for all s < — |. The value s = — | belongs to this 
interval. Finally, by the Sobolev embedding and Poincare inequality [2] , there are constants C, C > 
such that 

||«|li4 < C\\{u^y\\l2 < 4C7||n||ioo||n'||i2 < C\\n\\jj, < C\\u\\^^. 
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Since the norm in 'Hi(M) for the function u{x) is equivalent to the norm in l^^i"^) for the vector a by 
Lemma [H the estimate for the nonhnear vector field is completed by 

||g(a)||2 < Co||n||i4 < Collallf, , (2.9) 

-1/2 1/2 

for some Cq,Co > 0. The interpolation argument for the bounds (j2.7p and (|2.9|) concludes the proof 
that the nonlinear vector field f(a) + ag{a) maps ly2(^) to -^^2(N)- ^ 

Theorem 2 The discrete dynamical system /i2. 5\) is globally well-posed in the phase space a G /^y2(^)- 

Proof. By Proposition 2.2 in [6], the GP equation (jl.2p is globally well-posed in the phase space 
u G 'Hi(IR). By Lemma dl the trajectory u{t) G 'Hi(M) is equivalent to the trajectory a(i) G Z^^gC^^) 
on t G M. By Lemma [21 the vector field of the discrete dynamical system (|2.5p is well-defined on 
^1/2 (^) ^^(^)> where it is equivalent to the vector field of the GP equation ()1.2p by virtue of standard 
orthogonal projections. □ 



3 Existence and stability of stationary solutions 

Stationary solutions of the dynamical system (|2.5p take the form a(t) = Ae~*'^*, where A is a time- 

oo 

independent vector and /i is a parameter of the solution. If A G /^^(N) and = ^ An4>n{x), then 

n=0 

4> G 7^1 (M) is a stationary solution of the GP equation (|1.2p . that is (j){x) satisfies the ODE (|1.3p . The 
vector A is found as a root of the infinite-dimensional cubic vector field F : Z^^gl^) x K ^^1/2(^)5 
where the n-th component of F(A,/i) is given by 

Fn — {l^ ~ fl^-^n ~ C" ^ ^ -f^n;ni,n2,n3^rai^n2^n3 = 0, V?T. = 0, 1, 2, ... (3.1) 

(ni,n2,n3) 

The Jacobian operator DaF(0,/x) is a diagonal matrix with entries ^ — n and it admits a one- 
dimensional kernel if ^ = no for any non-negative integer no- The corresponding eigenvector is e„(j, 
the unit vector in /^(N). By the local bifurcation theory [8], each eigenvector of L'AF(0,no) can be 
uniquely continued in a local neighborhood of the point A = G t^i^i^) and ^ = no G M. We are 
particularly interested in the second eigenvalue no = 1, which corresponds to the dark soliton (j){x) with 
a single zero (node) at x = 0. (Other bifurcations of stationary localized solutions (j){x) are considered 
in [3].) Details of this bifurcation are given in the following proposition. 

Proposition 1 Consider real-valued roots (A,/i) of the vector field F(A,//) such that A G /^y2(^)- 
There exists a unique family of solutions near fJ. = 1 parameterized by e such that 

< C2e\ (3.2) 

for some e-independent constants Ci,C2 > and sufficiently small e. Moreover, if a ^0, the solution 
(A,/x) is smooth with respect to e for sufficiently small e and ■^Q{A) ^ 0, where Q{A) = ||A||p. 

Proof. Both F(A,/x) and Dji^F{A, /j,) are continuous in a local neighborhood of A = G /^^gC^) 
and /i = 1 G M. At the point A = and fi = 1, the operator has a one-dimensional kernel with 



eei||,2 < Cie 

'l/2 
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the eigenvector ei G /^(N). By using the method of Lyapunov-Schmidt reductions [S], we set A 



ei + A 



and /i = 1 + /i, where A is an orthogonal complement of ei in P{N) such that Ai = 0. The 



orthogonal projection of equation (|3.ip to ei gives a bifurcation equation for jl 

;l,ni,n2^ni^n2 ~^ 



1^ 



as 



ni 



E 

("1,"2) 



E 

(ni,n2,n3) 



AAA 

^ I;7ii,n2,n3^ni ^n2 



where the index for (ni, 712,71.3) in the summation signs runs on the set {0,2,3,...}. Let P be an 
orthogonal projection from to the orthogonal complement of ei. Then the inverse of PL>aF(0, 1)P 
exists and is a bounded operator from Z^gC^) to Z^gC^)- the Implicit Function Theorem, there 



exists a unique smooth solution A in the neighborhood of A = G Z^(N) such that ||A||^2 < Ci£ 

1/2 

for some Ci > 0. By the Implicit Function Theorem, there exists a unique smooth solution /i of the 
bifurcation equation in the neighborhood of /i = such that \fl — e^(T-fCi,i,i,i| < C2£^ for some C2 > 0. 



The value -RTi 1 1 1 



/32n 



is computed in Table I. Since Q(A) 



I Al|2 



1 



+ 0(e'^), then |;(5(A) / near = 1 for a / 0. 



+ 0(e4) and 
□ 





n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = b 


K 


1 


3 


41 


147 


8649 


323U7 


4n/2^ 


64 1/2^ 


256v^ 


16384 


65536v^ 


-f^l,n,n,l 








11 


75 


133 


2v^ 




16v^ 


32^2^ 


256v^ 


512v^ 




1 





1 

8v^ 





3v^ 





2v^ 


32 



Table I: Numerical values for Kn,n,n,n = \\(t>n\\li, -f^i,n,n,i = {4>iAn)i and ^^0,1,1," = ('/'o'/'n, </'i)- 

Let (A, /i) be a real- valued root of the nonlinear vector field (j3.ip such that A G /^^2(^)- Spectral 
stability of the stationary solution is studied with the expansion 

a(t) = e-*'^* [a + (B - C) e^^* + (B + C) e"^^* + 0(||B||2 + ||Cf)l , (3.3) 

where the spectral parameter G C and the eigenvector (B,C) G /^(N, C^) satisfy the linear problem 

L+B = l^C, L_C = 17B, (3.4) 
associated with matrix operators L±. Their ?2-th components are defined in the form 

(L+B)„ = {n- ll)Bn + 3fJ Yjni ^n,ni , 



(L_C) 



(n - fi)Cn + O-^ni K,niC'ni, 



Vn = 0,1,2,3,..., 



(3.5) 



where Ki 



"i,n2,n3^n2^n3 • have used here the symmetry of the coefficients Kn,ni,n2,n3 

{'n.2,n'i) 

with respect to the interchange of (711,712,71-3). 



ors 



Lemma 3 Let {A., /J-) be a real-valued root of the vector field F(A,^) such that A G ly^i^)- Operat 
L+ and L_ admit closed self-adjoint extensions in P{N) with the domain in Z^(N). 



Proof. The diagonal unbounded part of L± maps ifiN) to /^(N). We need to show that the non-diagonal 
part of L± represents a bounded perturbation from /^(N) to /^(N) if A G Z^^gl^)- This is done by using 
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the same ideas as in the proof of Lemma [2J 



E 

n=0 



ni 



n=0 (ni,n2,n3) (mi,m,2,m,3) 



Kn,n\ ,712 ,n3 Kn^rn\ ,m2 ,r(i3 ^n2 ^m2 -^ni -^mi 



U/2„||2 < |U,||4 
I " '^IIl2 ^ II "IIlo 



Il2 



n=0 



n=0 



where = I]^=o ^'^'^ ^(^) = I]jr=o ■^"■0n(x). 



□ 



Remark 2 The result of Lemma [3] is obvious from the equivalence between the space 7i2 (K) for the 
function v{x) and the space /f (N) for the vector B, see Remark [TJ We recall that the differential 
operators C± given by ()1.5p are defined on the domain 7i2{^) and the matrix operators L± given by 
(j3.5p represent the action of differential operators on the basis of Hermite functions in 7^2 ( 



The linear problem (j3.4p has eigenvalue i7 = of geometric multiplicity one and algebraic multiplicity 
two due to the exact solution 

L_A = 0, L+d^A = A, (3.6) 

where the smoothness of A with respect to /i near /i = 1 is guaranteed by Proposition [TJ 

When A = and ^ = 1, the spectrum of the eigenvalue problem (j3.4p is known in the explicit form. 
It consists of eigenvalues = and = ±1 of geometric and algebraic multiplicities two and simple 
eigenvalues = itm for all m = 2, 3, .... The double zero eigenvalue persists for any e according to the 
exact solution ()3.6p . stemming from the underlying U{1) invariance of the system. Splitting of all other 
eigenvalues in a local neighborhood of A = and /U = 1 is described by the following proposition. 



Proposition 2 Let (A,/x) he defined by Proposition [1\ for sufficiently small e. Non-zero eigenvalues 
of the linear problem ^3.4^ form a set {ibOm.}m=o simple real symmetric eigenvalue pairs, such that 



1| < Coe\ 



1 + 



8V27r 



and 



m + e'^a {Ki^i^i^i - 2Km+i,i,i,m+i)\ <Cme'^, Vm = 2,3, , 



(3.7) 



(3.8) 



for some e-independent constants Co,Ci,Cm > 0. 



Proof. Since the essential spectrum of the matrix operators L± is empty and the potential terms are 
bounded perturbations to the unbounded diagonal terms, isolated eigenvalues split according to the 
regular perturbation theory [10]. The formal power series expansion for a simple eigenvalue = m = 
2, 3, ... is defined by 

" B = e„+i +e2B + 0(e4), 
C = em+i+e^C + 0{e^), 
n = m + e'^n + Oie^). 

Projections to the component n = m + I lead to a linear system at the leading order O(e^) 



(3.9) 




cr [^1,1,1,1 
cr [^1,1,1,1 



3-f^m+l,l,l,m+l] + ^ 
-fi^m+l,l,l,m+l] + ^- 



(3.10) 
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The linear system has a solution if and only if Cl = cr (2A'm+i,i,i,m.+i — -^1,1,1,1) • Persistence of the 
eigenvalue by the perturbation theory results in the expansion (|3.8p . The power series expansion for 
the double eigenvalue = 1 is defined by 



= aeo+/5e2 + e2B + 0(e4), 
= -aeo + /3e2 + e2c + 0(e^), 
= l + e^n + 0{e^), 

where (a, /5) are arbitrary parameters. Projections to the components n = and n 
linear system at the leading order O(e^) 

^0 + Co) = (T [3Ko;i,i,oa + 3^0,1,1,2/3 - -fs^i, 1,1,1a] + Q,a 

' Co + Bo) = cr [i^o,i,i,oa - ^0,1,1,2/3 - 1,1,1a] + 

B2 -C2] = a [i^i,i,i,i/3 - 3i^2,i,i,oa - 3i^2,i,i,2/3] + 

'C2-B2) = a [i^i,i,i,i/3 + i^2,i,i,oa - ^2,1,1,2/3] + 

The linear system has a solution if and only if (a, (3) satisfies a homogeneous system 

a (i^i, 1,1,1a - 2Ko,i,i,oa - ^0,1,1,2/3) = (^a, 
(-i^i,i,i,i/3 + i^2,i,i,oa + 2i^2,i,i,2/3) = ^P- 



(3.11) 

2 leads to a 



(3.12) 



(3.13) 



The homogeneous system for (a, (3) has a non-zero solution if and only if satisfies a quadratic equation, 
roots of which are given by 



O = cr ( -fir2,i,i,2 - ^0,1,1,0 ± \ 



0,1,1,0 



K2,l,l,2Y - ^0,1,1,2 



(3.14) 



It follows from Table I that J (i^i, 1,1,1 — -f^o,i,i,o — ^2,1,1,2)^ — 1 2 



1 



1 

16v^' 



and -fir2,i,i,2--?^o,i,i,o = 
Persistence of the eigenvalues by the perturbation theory results in the expansion (j3.7p . □ 



16V27r 



Corollary 1 Let [Bm,Cm] he an eigenvector of the linear problem ^3.4\ ) for the eigenvalue £ 1^+ 
for any m = 0,1,2,3... in Proposition For sufficiently small e, the eigenvalue has positive 
signature o/(Bo,L+Bo), the eigenvalue ili has negative signature o/(Bi,L+Bi), while all other eigen- 
values rim with m = 2,3,... have positive signature of (Bm,-^^+Bm), where (•,•) denotes a standard 
inner product in P{H). 



Proof. In the case = 0, the homogeneous system (I3.13P for (a, (3) has a one-parameter family of 
solutions with (3 = — \/2a, such that (Bo,L+Bo) = — |ap + |/3p + O(e^) > for sufficiently small e. 
In the case 17 7^ 0, the homogeneous system ()3.13p for (a, (3) has a one-parameter family of solutions 
with a = — \/2/3, such that (Bi,L+Bi) = — |ap + + O(e^) < for sufficiently small e. In the case 
of other eigenvalues, it is obvious from the proof of Proposition [2] that (B^, L+B^) = m + O(e^) for 
m = 2,3,.... □ 

Remark 3 The double zero eigenvalue is associated with the expansions 



A = eei + 0(e3), = ^^ei + 0(e), (3.15) 

where e is sufficiently small. As a result, (A, d^A) = + O(e^). 
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Lemma 4 Let e be sufficiently small. Let {[B^, Cm]"^}m=o ^ of real-valued eigenvectors of the 
linear problem {3.4\) for the set of positive eigenvalues {^m}m=Q- eigenvectors is symplecti- 

cally orthogonal such that 

(B^.,C^) = 0, Vm'/m (B„,C„)/0, Vm = 0, 1, 2, 3... (3.16) 

In addition, two eigenvectors {[0, A]-^, [9^A, 0]-^} for the double zero eigenvalue 17 = are symplecti- 
cally orthogonal to other eigenvectors and {A,d^A) ^ 0. The set of eigenvectors 

{[Bm, C^f}^=o © {[B™, -Crnf}^=o © {[0, A]^, [d,A, of} (3.17) 
is a basis in /^(N, M^) which is orthogonal with respect to the symplectic projections 113. 16\) . 

Proof. All eigenvalues {Qm}m=Q positive and simple for sufficiently small e by Proposition [2l Since 
L± are self-adjoint in and 0^ is a real eigenvalue, then the eigenvector [Bm,Cm]"^ of the linear 

problem (j3.4p can be chosen to be real- valued. The orthogonality relations (I3.16P follow by direct 
computations from the linear problem ()3.4p for distinct eigenvalues ^l^n' for all m' ^ m. Values 

of (B,„,Cm) are proportional to the values of {Bjn,L+'Bm) for 0,^ they are non-zero for 
sufficiently small e by Corollary[TJ The value of (A, 5^A) = 5^Q(A) is non-zero for sufficiently small 
e by Proposition [TJ By the proof of Proposition [2] and Remark [3l the eigenvectors of the set (|3.17p are 
represented for sufficiently small e by the standard basis {em}m=o ® {^rn}m=o perturbed by a bounded 
perturbation in P(N) of the order O(e^). Also 

nm = m + 0{e'), {Bm, Cm) = ^^"^'^+^"^) = 1 + 0{m~'e^), Vm = 2,3,..., (3.18) 

for sufficiently small e, uniformly in m. Since no other eigenvalues exist, the set of linearly independent 
eigenvectors ()3.17p is complete in /^(N,M^). According to the Banach Theorem for non-self-adjoint 
operators, the set is a basis if and only if the spectral projections are bounded from below by a non- 
zero constant in the limit m ^ oo, which follows from the uniform asymptotic distribution (j3.18p . 
Therefore, the set (l3Tr|l is a basis in Z2(N,R2). □ 

Lemma 5 Fix e ^ sufficiently small. Simple positive eigenvalues of the set {Qm}m=o 
intervals 

a > : 1^0 = 1 and m - C„e^ <Qrn. <rn, Mm G N (3.19) 

and 

< : 9.0 = 1 and m<^m<m + C+e^^ Vm G N (3.20) 
for some e-independent constants > 0. 

Proof. The eigenvalue Qq = 1 persists for any e G M due to equivalence of the linear eigenvalue 
problems (jl.lip and (|3.4p for G Hi{R) and A G ll/2(^) ^^'^ existence of the exact solution (|1.12p 
of the linear eigenvalue problem (jl.lip . The corresponding eigenvector (B, C) of the linear eigenvalue 
problem (j3.4p is found from the series representation 

oo oo 
= ^ BnCpnix), -X(/)(x) = ^ C„(/>„(2;). 

n=0 n=0 
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The eigenvalue ili satisfies the bounds (|3.19|) - (j3.20|) due to the explicit bound (|3.7|) . We use the 
bound (jS.Sp to prove the bounds (j3.19p - ()3.20p for eigenvalues for all m = 2,3,.... The values of 
-^1,1,1,1 ~ 2i^m_(_i 1 1 m+i are positive for the first values of m = 3,4, ... as follows from Table I, e.g. 

21 _ 59 



-f^l, 1,1,1 — 2^3,1,1,3 = , -fi^l, 1,1,1 - 2i^4,i,i,4 = ———=, i^l, 1,1,1 — 2X5,1,1,5 



16V27r 128V27r 256V27r 



Note that the positive numerical values are monotonically increasing. According to the main theorem 
in [7], the sequence {||(/'n||L4}„gN is monotonically decreasing to zero with the bound (|2.8p . Since 
i^m+i,i,i,m+i < Il0i|li4||0m+i|li4, then 



-f^l, 1,1,1 - 2-R'm+l,l,l,m+l > 110111^4 (||0l|li4 - 2||(/)m+i|||4) , Vm = 2,3, ... 

Since ||(/)m+i||^4 decays monotonically to zero as m ^ oo, there exists M sufficiently large, such that 
the lower bound above is strictly positive for m > M. □ 



4 Existence of periodic solutions 

Let (A,/i) be a real-valued root of the nonlinear vector field F(A,/i) such that A G /^^gl^)- We use a 
decomposition a(t) = e~*^* [A + B(t) + iC{t)] with real-valued vectors B and C to rewrite the discrete 
dynamical system (|2.5p in the form 

B = L_C + (tN_(B,C), -C = L+B + (jN+(B,C), (4.1) 

where the operators L± are defined by (j3.5p and the vector fields N-|-(B,C) contains quadratic and 
cubic terms with respect to (B, C). By Theorem [2l the initial- value problem for system (j4.ip is globally 
well-posed and the solution set (B, C) G /|^2(N)I^^) is equivalent to the solution set {v,w) G 'Hi(R,M^) 
of the PDE system (jl.4p . The discrete dynamical system ()4.ip inherits the Hamiltonian function (jl.9p 
in the form 

1 1 

H = -(B,L+B) + -(C,L_C)+(J Yl ^" 

(n,ni,n2,n3) 

~l~ ~^ ^ ^ Kn^n\,n2,n^ iBn^Bn2Bfi-^BYi + 1B^-^Bji2Cn-:^Cn -\- Cj^j^ Cjjj C'ng C^) (4'2) 
(n, 711,712, 713) 

and the conserved quantity (jl.lOp in the form 

Q = 2(A, B) + (B, B) + (C, C). (4.3) 

Using LemmaSl we represent a solution (B, C) of the discrete system ()4.ip by the series of eigenvectors 
()3.17p associated with the linear problem (|3.4p : 



B(t) = E^=o&-WBm + E^=o^'-(*)Bm + /5(t)a^A, 



where 6o(0) t)(t) = (6i,625---) are complex-valued and /9(t), 7(t) are real-valued. The linear part of 
system (14. ip becomes block-diagonal in the representation (14. 4|) . yielding the evolution equations 

6m - ■i^^TTi^m = o-A'^m(&o,b,/3,7), Vm = 0,1,2,3... (4.5) 
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and 

/3 = aSo(6o,b,/3,7), j + P = aSi{bo,h, (4.6) 

where 

^\*-^m) -Dm/ 

and 

(A,N_(B,C)) ^ , ^ (a^A,N+(B,C)) 
5o(6o,b,/3,7) = ^A,0,A) ' ^i(^o,b,/3,7) = (a,0,A) " 

Using conservation of Q given by (|4.3p and the decomposition ()4.4p . one can integrate the first equation 
of system (I4.6P in the form 

where Q is constant in time t £ M. As a result, the second equation of system (j4.6p is rewritten 
exphcitly in the form 

. _ ||B||2 +||C||2, -2a(9^A,N+(B,C))-g 

We are now ready to apply the method of Lyapunov-Schmidt reductions to the proof of Theorem [TJ 

Proof of TheoremUl The vector space (B, C) S /^y2(^'^^^) is equivalent to the vector space b G /^^gC^) 
because of the asymptotic distribution (|3.18p . For instance, one obtains that 

oo oo 

+ n)\Bn\ ~ (B, L+B) = 2 ^m{Cn„ B^)|6^p + \(3\^{A, d^A) ~ + 

n=0 m=0 neN 

We should work in the space of T-periodic functions 6o(t), h{t) G ly^i^)^ Pi't) and j{t) on t G M, 
where T is close to 27r. This period corresponds to the eigenvalue = 1 which persists for any e G M. 
By Lemma EJ all other eigenvalues of the linear problem (j3.4l) satisfy the non-resonance conditions 
n 7^ Urn, Vn, m G N for any fixed e ^ sufficiently small. As a result, we define periodic functions b(t), 
and ^{t) in terms of the periodic function bQ{t), which solves a reduced evolution problem. Let 6 
be sufficiently small. We shall prove that there exist solutions of system (14. Sp . (14. 7p and (14. 8p which 
are T-periodic on t G M satisfying the apriori bounds 

\bo{t)\ < e6Co, \\Ht)\\i2 ^ < e6^Cb, \P{t)\ < e^S^Cfs, h{t) - 6a\ < e'^S^C^, Vt G M, Va G M, (4.9) 

1/2 

for some (e, (5)-independent constants CQ,Ch,Cp,C^ > 0. If 6o(t), h{t) G Z^^2(^)i Pi't) and 7(t) are T- 
periodic functions on t G M satisfying the bounds diJ]), then (B(t), C(i)) G /^^2(^'^^) 

is a T-periodic 

function on t G M satisfying the bound 

l|B(t)||i2^^ + ||C(t)||;2^^ < C£6, Vi G M, Va G M, (4.10) 

for some (e, (5)-independent constant C > 0. Here we recall the expansion (j3.15p for A, 5^A and the 
fact that (Bm,Cm)"^ are close to the unit vectors for sufficiently small e. Since N-|-(B,C) is cubic 
with respect (A,B,C), contains quadratic terms in (B,C), and maps Z^^gC^i^^^) to '?.i/2(^' ^^)' 
obtain the bound 

||N±(B(t),C(t))||i2 <C±e^6'^, VtGM, VaGM, (4.11) 

— 1/2 
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for some (e, 5)-independent constants C± > 0. By the Implicit Function Theorem to the right-hand- 
side of equation (j4.8p . there exists a unique constant Q in the interval \Q\ < Cqe^S'^ for some Cq > 0, 
such that the periodic function in the right-hand-side of equation (j4.8p has zero mean on t G M. In this 
case, there exists a periodic solution 7(t) = 6a + j{t) of the differential equation (j4.8p . where ^{t) is 
a uniquely defined varying part and 6a is an arbitrary mean part. The varying part 7(t) satisfies the 
last bound in the list (j4.9p . The function P{t) is uniquely defined by the explicit representation ()4.7p 
and it hence satisfies the third bound in the list (j4.9p . 

Consider now system (14. Sp for m G N. Recall that 0^ —rn = O(e^) for m = 1, 2, ... uniformly in m G N 
for sufficiently small e. By the Implicit Function Theorem, there exists a unique solution h{t) G ll/^i^) 
defined by the periodic function 6o(^) and parameter a G M for sufficiently small 6 provided that the 
distance \^rn — 7^ and the frequency of the periodic function 60 (t) is such that 17 — > 1 as (5 — > 0. 
By the bound (j4.1ip and the distribution Vim — m = O(e^) for all m G N, the function h(t) satisfies the 
second bound in the list ()4.9p . 

Eliminating the components b, /5 and 7 from equation ()4.5p for n = 0, we obtain a reduced evolution 
problem for 60 (t) in the form 

60 = ^60 + ^(^0; a), (4.12) 
where R{bo;a) is a remainder term. Explicit computations of NQ{bo,h, (3,^) show that 

R{bo;a) = e[iKi{£)bl + iK2i£)bl + iK3{e)\bo\'^ + iK^{e)6'^a'^ + K5{e)6abo] 

+ 0{\bo\\eH^a^\bo\,e\bo\\\h\\), (4.13) 

where i^i,2,3,4,5 are real-valued constants which are bounded for sufficiently small e. We are looking 
for T-periodic functions 60 (t) which satisfy the evolution problem (j4.12p . have the leading order 60 ~ 
Sf^e'*"'"*'^ , where r G M is arbitrary, and satisfy the first bound in the list ()4.9p . By the normal form 
analysis of the ODE (|4.12p (see [H]), the quadratic terms in the remainder ()4.13p do not change the 
frequency of oscillations of the periodic function 60 (t) at the leading order and therefore, — 1| < 
Cqe^iJ^ for some Cq > 0. Since the Hamiltonian function (j4.2p of system (j4.ip is constant in time, it 
remains constant when the function 60 (t) solves the reduced evolution problem (I4.12P and the functions 
b(i), (3{t) and 7(t) are constructed above. By the normal form analysis of reversible systems, there 
exists a two-dimensional invariant manifold of system ()4.12p filled with periodic solutions of frequencies 
close to Q = 1 and parameterized by {6,t) in addition to parameter (e,a). □ 

Remark 4 Theorem [1] is reminiscent of an infinite-dimensional analogue of the Lyapunov Theorem 
for persistence of periodic orbits in Hamiltonian systems (see Chapter II, Section 45 on pp. 166-180 
of [E]). However, due to the symmetries, a double zero eigenvalue occurs in the linear problem (j3.4p . 
and the proof of Theorem [1] is complicated by the analysis of the associated two-dimensional subspace. 
Similar theorems on persistence of fc-dimensional tori in n-dimensional Hamiltonian system with k — 1 
additional conserved quantities were studied in the Nekhoroshev-Kuksin Theorems (see Theorem 2.3 
on p. 4 of [1] and Theorem 1 on p. xiii of [13]). 

Remark 5 The periodic solution of Theorem[T]has the smallest frequency in the focusing case a = — 1, 
since J^i > 1 in the bound (j3.20p for sufficiently small e. However, it is not the smallest frequency in 
the defocusing case a = 1 since J7i < 1 in the bound (j3.19p . Persistence of the periodic solution for 
the smallest frequency f^i can not be proved by a simple application of the Lyapunov Theorem since 
the bound (|3.19p does not guarantee that the non-resonance conditions nJ7i 7^ ai'e satisfied for all 
n G N and m = 2,3, .... By the same reason, persistence of quasi-periodic oscillations on the tori with 
two and more frequencies {l,^}i,0,2, ■■■} can not be proved for small e. 
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Remark 6 Persistence of quasi-periodic oscillations on the tori along the Cantor set of parameter 
values was proved in Section 2.5 on p. 33 of [l3_j for the Hartree nonlinear functions and a perturbation 
of the parabolic potential V{x) = ^x^ by a localized potential Vo{x). Our main result is stronger 
than this application of the main theorem in [13] since the periodic orbit is continuous with respect to 
parameters of the PDE problem rather than along the Cantor set of parameter values. 

5 Numerical Results 

We illustrate results of our manuscript with some numerical approximations. First, we identify the 
relevant branch of stationary solutions of the ODE (jl.3p . To do so, we use a fixed point method 
(Newton-Raphson iteration) to solve a discretized boundary-value problem. A centered-difference 
scheme is applied to the second-order derivatives with a typical spacing Ax G [0.025,0.1]. We are 
using a sufficiently large computational domain x G [— L, L] such that the boundary conditions do not 
affect the approximations within the considered numerical precision. The solutions (p{x) are obtained, 
using continuation, as a function of parameter /x. The continuation of the solution branches is per- 
formed from the linear limit /i = 1, both for the cases a = 1 and cr = — 1. The results are shown 
in Figure [H illustrating the quantity Q = \\<j)\\'^2 as a function of //. The numerical findings are also 
compared to the asymptotic result ()3.2p of Proposition [T] indicating the good agreement of the latter 
prediction with our computational results for a fairly wide parametric window. 

Once the corresponding numerical solution is identified (for a given a and fi), the linear eigenvalue 
problem (jl.lip is approximated numerically. We use again a discretization of differential operators on 
a finite grid, such that the spectral problem (jl.lip becomes a matrix eigenvalue problem that is solved 
through standard numerical linear algebra routines. The relevant lowest eigenvalues are presented in 
Figure [2] and are also compared with the corresponding asymptotic results (|3.7p - (j3.8p of Proposition 
[2j The dashed lines show asymptotic results (IA.7p - ()A.8p of Appendix A derived in the limit fi ^ oo 
for a = 1. Once again, the good agreement offers us a quantitative handle on the relevant eigenvalues. 

Finally, we have also examined periodic oscillations of dark solitons in the numerical simulations of 
the GP equation (jl.2p . A typical example is shown in Figure [3] for a = 1 and /i = 1.1 for the initial 
condition n(x,0) = (j){x) + 6(j)'{x) with 5 = 10^^. The top left panel shows the space-time contour 
plot of |ii(x,f)p, clearly highlighting that this is a small (imperceptible, at the scale of this panel) 
perturbation of a stable stationary solution (j){x). The bottom left panel shows the space-time contour 
plot of \u{x, t)p — (/>'^(x), emphasizing the time-periodic oscillations of the perturbation to the stationary 
solution. The periodic oscillations are also visible on the top right panel where |n(xo,i)P is plotted 
versus t for xq = 2. Finally, the bottom right panel illustrates the Fourier transform of the time series 
of |ti(xo,t)P (normalized to its maximum). It shows a high peak of the frequency spectrum near the 
value = 1 , in agreement with the results of the main Theorem [TJ 

Acknowledgement. D.P. thanks to W. Craig and V. Konotop for useful discussions related to the 
project. D.P. is supported by the Humboldt and EPSRC fellowships. P.G.K. is supported by NSF 
through the grants DMS-0204585, DMS-CAREER, DMS-0505663 and DMS-0619492. 

A Asymptotic distribution of eigenvalues 

Let us consider the case a = 1, when the solution </)(x) of the ODE (jl.3p bifurcates to the interval 
^ > 1 (see Proposition [T] and Figure [T]). We are interested in the distribution of eigenvalues of the 
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Figure 1: Branches of dark solitons versus fx both for the case of o" = —1 (when fj, < 1) and a = 1 
(when ^ > 1). The numerically obtained solution is shown by solid line and the asymptotic solution 
(|3.2p is shown by dash-dotted line. 




^^0 ^5 |i 5 10 15 



Figure 2: Smallest purely imaginary eigenvalues of the linear eigenvalue problem (11. lip versus ^. The 
numerically obtained eigenvalues are shown by solid lines, the asymptotic results (|3.7p - (|3.8p are shown 
by dash-dotted lines, and the asymptotic results (|A.7p - ()A.8p are shown by dashed lines. 
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Figure 3: A typical example of the robust time-periodic solution of the Gross-Pitaevskii equation ()1.2p 
for cj = 1, /X = 1.1 and n(x,0) = (p{x) + 6<p'{x) with 5 = 10~^. The top left panel shows the space-time 
contour plot of |n(x,t)p, the bottom left panel shows the space-time contour of |ti(x,t)p — <fP'{x). The 
top right panel shows the time evolution of |n(xo,t)P with xq = 2, while the bottom right panel shows 
the Fourier transform of the time series of |n(xo,t)P, featuring a peak at ~ 1. 
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linear problem p.ll|) as // — > oo, assuming that the solution 4>{x) persists in this hmit. It follows from 
the scaling transformation below equation (jl.2p that the limit /i — > oo of the normalized equation (jl.2p 
corresponds to the limit e ^ in the original GP equation We shall replace + ^ = /i and drop 

tilde notations for the sake of simplicity. We report here formal results based on asymptotic methods. 
Rigorous justification of these results is beyond the scope of our work. 

Denote the ground state of the ODE ()1.3p by (j)o{x) such that </>o(x) is even and positive on x G M 
and it decays to zero as |x| — > oo sufficiently fast. Using the substitution (j)o{x) = \J /ug(^) and ^ = -^^i 
we obtain an equation for 

which is solvable with the nonlinear WKB series [12]. The main result of the formal WKB theory is 
that there exists a classical solution q^i^C) of the ODE (jA.ip for sufficiently large /u > 1 such that 

= ^|«|;; (A.2) 

The linear problem (jl.lip associated with the ground state <^q(x) for o" = 1 and /i + | — > ^ can be 
written in variables v = V{^), w = W{(,) and A = for sufficiently large > 1. In new variables, it 
takes the form 

L+V = -AW, L^W = AV, (A.3) 

where 

L, = 3<,(e)-l+«=-jij-|, (A,4) 

Eliminating V{x), we close the linear problem (|A.3P at the fourth-order ODE 

L+L^w = rw, r = -A^. (A.5) 
By using the WKB theory ()A.2p . we consider the auxiliary eigenvalue problem 

1 p^{iv) _ i^-p ^r" = rw{0, ve e [-1, 1], (A.6) 



16//4 2/i2 

for W £ 1,1]). The entire spectrum of the problem ()A.6P is defined by a set of polynomial 

solutions W = Pm{€) = C™" + am,m-2'^"'~^ + ••• + Oim,ki^-, Vm e N, where A; = 1 if m is odd and /c = 
if m is even. The balance of the largest term in the ODE (|A.6p shows that the eigenvalue T = F.^ is 
found explicitly as T^ = "^''^^^ , while all coefficients {am,m-2k}^k^i^ are uniquely defined. Converting 
the values of F to the values of A, we have found that the linear problem (jl.lip associated with the 
ground state (j)o{x) has a set of simple purely imaginary and symmetric eigenvalue pairs {iirij^jInisNj 
such that 

lim = g (A.7) 

M^oo V2 

in addition to the double zero eigenvalue A = 0. 

Finally, the dark soliton (p{x) of the ODE (jl.Sp is obtained asymptotically from the ground state 
(l)o{x) by the factorization (j){x) = (I)q{x)iP{x), where ijj{x) is odd on x G R, positive on x € and may 
approach to the constant values as \x\ — > oo |16j . Using this factorization and the formal asymptotic 
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analysis, it was shown in [.16J that the spectrum of the hnear problem associated with the dark 

soliton (l>{x) admits a pair of simple purely imaginary eigenvalues such that 



lim ^0 = ^. (A.8) 

Although the analysis of [TU] was directed to the original GP equation in the limit of small e 

and the eigenvalue pair was found to be A — > itie, the scaling transformation to the normalized GP 
equation ()1.2p implies that A = -^^^ 

Numerical computations (see Figure [2]) suggests that the entire spectrum of the linear problem (jl.lip 
associated with the dark soliton (j){x) is a superposition between an infinite set of eigenvalues (IA.7|) of 
the linear problem (jl.lip associated with the ground state <j)o{x) and the additional pair of eigenvalues 

(lAl). 

Note that the linear eigenmode corresponding to the smallest eigenvalue = ™ay not result 
in the periodic solution of the nonlinear PDE system ()1.4p because the non-resonance condition n ^ 
\J m{m + 1) for all n,m S N is violated in the limit n,m ^ oo. Similarly, the linear eigenmode 
corresponding to the second eigenvalue = 1 may not result in the periodic solution of the PDE 

system (jl.4p because the non-resonance condition n ^ V™^+^ ). £qj. g^^j n,m = 2,3,... is violated at 
least for n = 6 and m = 8. In both cases, the Lyapunov Theorem for persistence of periodic orbit in 
Hamiltonian dynamical systems can not be applied [15\ - 

References 

[1] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, 
and Mathematical Tables (Dover, New York, 1965), chapter 22. 



[2: 



[4 



[6; 



R.A. Adams, Sobolev Spaces (Academic Press Inc., San Diego, 1978) 

G.L. Alfimov and D.A. Zezyulin, "Nonlinear modes for the Gross-Pitaevskii equation - demon- 
strative computation approach", arXiv: ^nlin. PS/0703006^ 

D. Bambusi and G. Gaeta, "On persistence of invariant tori and a theorem by Nekhoroshev" , 
Math. Phys. Electr. Journal 8, paper I (2002) 

V.A. Brazhnyi and V.V. Konotop, "Evolution of a dark soliton in a parabolic potential: application 
to Bose-Einstein condensates". Physical Review A 68, 043613 (2003) 

R. Carles, "Remarks on nonlinear Schrodinger equations with harmonic potential", Annales Henri 
Poincare 3, 757-772 (2002) 

G. Freud and G. Nemeth, "On the Lp-norms of orthonormal Hermite functions", Studia Scien- 
tiarum Mathematicarum Hungarica 8, 399-404 (1973) 

M. Golubitsky and D.G. Schaeffer, Singularities and Groups in Bifurcation Theory, vol. 1, 
(Springer- Verlag, New York, 1985) 

B.L.G. Jonsson, J. Frohlich, S. Gustafson, and I.M. Sigal, "Long time motion of NLS solitary 
waves in a confining potential", Annales Henri Poincare 7, 621-660 (2006) 



19 



[10] T. Kato, Perturbation theory for linear operators, (Springer- Verlag, New York, 1976) 

[11] V.V. Konotop, "Dark solitons in Bose-Einstein condensates: theory" in "Emergent Nonlinear 
Phenomena in Bose-Einstein Condensates", Eds. P.G. Kevrekidis, D.J. Franzeskakis, and R. 
Carretero-Gonzalez (Springer- Verlag, New York, 2007) 

[12] V.V. Konotop and P.G. Kevrekidis, " Bohr-Sommerfeld quantization condition for the Gross- 
Pitaevskii equation". Physical Review Letters 91, 230402 (2003) 

[13] S.B. Kuksin, Nearly Integrable Infinite-Dimensional Hamiltonian Systems (Springer- Verlag, 
Berlin, 1993) 

[14] Y.A. Kuznetsov, Elements of Applied Bifurcation Theory, 2nd ed., Appl. Math. Sci. 112 (Springer- 
Verlag, New York, 1998) 

[15] M.A. Lyapunov, The General Problem of the Stability of Motion (Taylor and Francis, London, 
1992) 

[16] D.E. Pelinovsky, D. Frantzeskakis, and P.G. Kevrekidis, "Oscillations of dark solitons in trapped 
Bose-Einstein condensates", Physical Review E 72, 016615 (2005) 

[17] D.E. Pelinovsky and P.G. Kevrekidis, "Dark solitons in external potentials", Zeitschrift fiir Ange- 
wandte Mathematik und Physik, to be published (2007) 

[18] N.G. Parker, N.P. Proukakis, C.F. Barenghi, and C.S. Adams, "Dynamical instability of a dark 
soliton in a quasi-one-dimensional Bose-Einstein condensate perturbed by an optical lattice" , Jour- 
nal of Physics B: Atomic Molecular Optical Physics 37, S175-S185 (2004) 



20 



